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ABSTRACT 



We describe the HST ACS Coma cluster Treasury survey, a deep two- 
passband imaging survey of one of the nearest rich clusters of galaxies, the Coma 
cluster (Abell 1656). 

The survey was designed to cover an area of 740 arcmin^ in regions of differ- 
ent density of both galaxies and intergalactic medium within the cluster. The 
ACS failure of January 27th 2007 leaves the survey 28% complete, with 21 ACS 
pointings (230 arcmin^) complete, and partial data for a further 4 pointings (44 
arcmin^) . 

Predicted survey depth for lOcr detections for optimal photometry of point 
sources is g' = 27.6 in the F475W filter, and 1^=26.8 mag in F814 (AB magni- 
tudes). Initial simulations with artificially injected point sources show 90% re- 
covered at magnitude limits of g' = 27.55 and Ic = 26.65. For extended sources, 
the predicted lOa limits for a 1 arcsecond^ region are g' = 25.8 mag arcsec"^ and 
Iq = 25.0 mag arcsec"^. 

We highlight several motivating science goals of the survey, including study of 
the faint end of the cluster galaxy luminosity function, structural parameters of 
dwarf galaxies, stellar populations and their effect on colors and color gradients, 
evolution of morphological components in a dense environment, the nature of 
ultra compact dwarf galaxies, and globular cluster populations of cluster galaxies 
of a range of luminosities and types. This survey will also provide a local rich 
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cluster benchmark for various well known g'/o 6a/ scaling relations and explore new 
relations pertaining to the nuclear properties of galaxies. 

Subject headings: galaxies: clusters: individual (Abell 1656) — galaxies: fun- 
damental parameters — galaxies: photometry — galaxies: stellar content — 
galaxies: structure — galaxies: luminosity function 

1. Introduction and background 

The Coma cluster of galaxies, Abell 1656, is along with the Perseus cluster the nearest 
rich, and dense cluster environment. Unlike the Perseus cluster, it is at high galactic latitude 
(b = 87.9°) and it has been a popular target for study at all wavelengths. Progressively 
deeper wide-area photometric surveys of Coma have become available over the past 30 years, 
and waveband coverage has spread from the original B and V band surveys into the near 
ultra-violet and infra-red (Godwin & Peach 1977; Godwin, Metcalfe & Peach 1983 (GMP); 
Komiyama et al. 2002; Adami et al. 2006; Eisenhardt et al. 2007). A larger area around Coma 
is covered by the imaging part of Data Release 5 of the Sloan Digital Sky Survey (Adelman- 
McCarthy et al. 2007). From these surveys, samples of galaxies have been selected for 
spectroscopic study, which has resulted in an understanding of the internal dynamics of the 
cluster (CoUess & Dunn 1996; Mobasher et al. 2001; Edwards et al. 2002; Gutierrez et al. 
2004), the internal dynamics of cluster members (Lucey et al. 1991; Jorgensen et al. 1996; 
Smith et al. 2004; Matkovic & Guzman 2005; Cody et al. 2007), and their mean luminosity 
weighted stellar ages, abundances and a-enhancement (Bower et al. 1992a,b; Guzman et 
al. 1992; Caldwell et al. 1993; Jorgensen 1999; Poggianti et al. 2001; Moore et al. 2002; 
Nelan et al. 2005; Smith et al. 2006; Matkovic et al. 2007). Coma presents us with the 
best opportunity to study large samples of galaxies of different luminosity, environment and 
morphological type, but at a common distance, and with a common Galactic extinction. 

There is good agreement on the distance to the Coma cluster, with independent studies 
using six different methods yielding values in the range 84 - 108 Mpc, as summarised in 
Table [H These values fit well with the current concordance cosmology: assuming Hq = 71 
km/s/Mpc; Q\ = 0.73; Qm = 0.27, and a redshift z = 0.0231 gives a distance of 99.3 Mpc. 
In this paper we assume a distance of 100 Mpc, equivalent to a distance modulus of 35.00. 
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Table 1. The Distance to Coma determined by different techniques 



Technique 


Distance (Mpc) 


Distance Modulus 


Reference 


I-band TuUy-Fisher 


86.3±6 


34.68±0.15 


TuUy & Pierce (2000) 


K'-band SBF 


85±10 


34.64±0.27 


Jensen et al. (1999) 


I-band SBF 


102±14 


35.04±0.32 


Thomsen et al. (1997) 




96±6 


34.90±0.14 


Gregg (1997) 


Fundamental Plane 


108±12 


35.16±0.25 


Hjorth & Tanvir (1997) 


Globular Cluster LF 


102±6 


35.05±0.12 


Kavelaars et al. (2000) 
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Coma lies in a rich region of the large-scale distribution of galaxies, at the intersection 
of a number of filaments. Figure [1] shows two projections of the distribution of galaxies 
in supergalactic co-ordinates in the region of Coma and the nearby richness class 2 cluster 
Abell 1367. The Great Wall (Geller & Huchra 1989) a vertical structure in the two panels 
of Figure [H runs through these two clusters, other filaments intersect the Great Wall at the 
Coma cluster. 
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Fig. 1. — The location and environment of the Coma cluster in supergalactic co-ordinates. 
The left panel shows a projection onto supergalactic Y, close to the plane of the sky. The 
right panel shows a projection onto supergalactic X, in this panel the horizontal (Y) axis 
is close to measured redshift. All axes are in equivalent cz. Positions of the points are 
derived from measured sky positions and redshifts. Members of Coma (centre) and Abell 
1367 (bottom) are plotted as filled red circles, for these galaxies the redshift used to compute 
the distance is the cluster redshift, with a random velocity offset chosen to make the cluster 
appear round. For non cluster members, measured cz is used, and the positions are plotted 
as filled black circles if they have velocities within ±600 km/s of the Coma mean in the depth 
dimension in that particular panel, and as open grey circles otherwise. The Great Wall is 
the structure at SGY ^ 6800 km/s in the right panel and is seen face on in the left panel 
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There is a uniquely rich multi-wavelength dataset on the Coma cluster. X-ray observa- 
tions covering a large area of the cluster have been made with ROSAT (White et al. 1993) 
and XMM-Newton (Briel et al. 2001), showing the distribution and properties of the hot 
intra-cluster medium (ICM), and X-ray properties of individual galaxies have been studied 
by Finoguenov ct al. (2004) and by Hornschemeier ct al. (2006). Coma has been shown 
by INTEGRAL to be an extended hard X-ray/soft 7-ray source (Renaud et al. 2006). At 
soft X-ray and Extreme Ultraviolet wavelengths there is a thermal excess (Lieu et al. 1996; 
Kaastra et al. 2003; Bonamente et al. 2003; Bowyer et al. 2004). GALEX has been used 
to observe the cluster in the mid and near Ultraviolet and has sufficient spatial resolution 
to measure the UV properties of individual galaxies. In the infra-red, studies of the galaxies 
have been made using SPITZER, both with IRAC at 3.6 - 8 //m (Jenkins et al. 2007) and 
with MIPS at 24 and 70 /im (Bai et al. 2006). At radio wavelengths, VLA continuum maps 
cover much of the cluster (Miller et al. in preparation), and in the HI 21 cm line there are 
extensive VLA imaging surveys (Bravo- Alfaro ct al. 2000, 2001), and single-dish spectra and 
fluxes for samples of spiral galaxies (Gavazzi et al. 2006; Vogt ct al. 2004). 

This wealth of existing data makes Coma a prime target for studies of the origin and 
evolution of the galaxy content of clusters, and of its interaction with the other components 
(gas and dark matter). Moreover, as the richest and best studied local cluster. Coma is the 
zero-redshift baseline for many studies of high-redshift clusters (e.g. Jorgensen et al. 2006). 
Comparison between low- and high-redshift clusters is vital for our understanding of their 
evolution, which in turn is essential if we are to disentangle evolutionary effects from the 
properties which tell us about their formation. 

We describe here the HST/ACS Coma Cluster Treasury survey, which aims to provide an 
unparalleled database of high spatial resolution images of a sample of cluster galaxies. At the 
distance of the Coma cluster (~100 Mpc), the resolution of HST/ACS (0".l) corresponds to 
~50 pc. This gives essentially the same physical resolution as ground-based observations have 
in Virgo and Fornax. Thus the HST/ACS Coma database provides a valuable comparison 
between high- and low-density clusters, for studies of the effect of environment on galaxy 
components. 

Whilst the HST observations are the prime data upon which this survey is based, it 
has already generated numerous ancillary observations with facilities such as Subaru, Keck, 
MMT, UKIRT. It is concurrent with surveys of the cluster in other wavebands, including 
the ultra-violet (GALEX), infra-red (SPITZER), X-ray (Chandra and XMM-Newton) and 
radio (VLA). 
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2. Scientific Motivation 

2.1. The Luminosity Function 

The logarithmic slope of the low-mass end of the Cold Dark Matter (ACDM) mass 
function has a slope a ^ —1.8 (e.g. Somerville & Primack 1999). In contrast, the faint end 
of the field-galaxy luminosity function has a slope a ~ —1.3, when measured either from 
optically-selected surveys (Blanton ct al. 2005), or Hi (Zwaan ct al. 2005) or a ~ —1.0, when 
measured in the K-band (Gardner ct al. 1997; Kochanck ct al. 2001). Luminosity functions 
in clusters and groups are often not well fit by a single Schechter function (Smith et al. 1997), 
and arc better modeled by a combination of a Gaussian and a Schechter function (Ferguson 
& Sandage 1991). The composite behavior of the LF and the trends with cluster richness are 
beginning to be understood in the context of the conditional luminosity function (Cooray 
& Cen 2005), which provides a powerful conceptual framework for exploring the physics of 
galaxy formation via studies of halos (clusters and groups) of different mass. The suppression 
of the faint-end slope of the LF relative to the CDM mass function is widely believed to be 
due to photoionization by the meta-galactic UV background, which suppresses star-formation 
in low-mass halos. This predicts an environmental dependence (Tully et al. 2002) because 
the fraction of dwarf-mass subhalos that collapse before re-ionization is larger in higher mass 
halos, i.e., much larger in a IO^^Mq halo than in a IO^^Mq halo. This prediction can be 
tested directly with our survey, as the very faint-end of the LF slope in Coma should be 
closer to the CDM prediction, and at least some of the dwarfs should have very old stellar 
populations. 

From earlier HST/WFPC2 data of a small area around NGC 4874 in the core of the 
Coma cluster, Milne et al. (2007) determine a steep faint-end slope (a ~ —2.0), which agrees 
with the slope for the faintest objects in a small-area imaging study by Bernstein et al. 
(1995), but not with studies using spectroscopic redshifts (Mobasher et al. 2003). 

Measurements of the faint-end LF slope in the Virgo cluster have been made by a 
number of authors, e.g. Trentham & Hodgkin (2002), and Sabatini et al. (2007). Trentham 
& Hodgkin (2002) find that the LF rises rapidly {a ^ —1.6) between Mb — —16 and 
Mb — —14, but fiattens faintward of this. The results of Sabatini et al. (2007) again suggest 
a steep slope between Mb = —16 and Mb = —14, fiatter between Mb = —14 and Mb = —11. 
Trentham & Tully 2002 analyse the slope of the faint end of the LF in a number of lower 
density groups. They find a mean slope of o; ~ —1.2 over a large range from Mb = —18 to 
Mb = —10. Micskc ct al. (2007) use SBF measurements to determine cluster membership in 
Fornax, and find an even fiatter faint-end slope, similar to that in the local group {a ~ —1.1). 
These authors point out however that the absolute number of dwarfs per unit virial mass is 
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higher in the lower density environments. Sabatini et al. (2007) concur that the luminosity 
function slope in the field is flatter than in the Virgo cluster, although they parameterise 
this in terms of a Dwarf to Giant Ratio (DGR). These studies present a consistent picture 
of a faint-end LF slope which is strongly dependent upon the density of the environment, 
and our study of the Coma LF, at a variety of galaxy densities, will present an important 
extension to this work. 

The area of the survey reduces the vulnerability of determinations of the faint-end slope 
to poor statistics and cosmic variance (Driver ct al. 1998), and provides a sufficient range of 
galaxy density to test variations in the faint-end LF across different cluster environments (e.g. 
Phillipps et al. 1998). The depth and spatial resolution allow the use of surface brightness 
and morphological criteria to assess the probability of cluster membership (Trentham et al. 
2001; Trentham & TuUy 2002) and examination of the relative contribution of dE, dSO and 
dirr galaxies to the faint end LF, and of the possible evolutionary relationship between these 
types (Aguerri et al. 2005; Mastropietro et al. 2005). 

Different formation mechanisms give rise to the disk, bulge and bar components of 
galaxies. To understand the relative dominance of these physical processes in the Universe 
therefore requires one to construct not simply a galaxy luminosity function, but separate 
bulge, bar and disk luminosity functions (Driver et al. 2007a; Laurikainen et al. 2005). 
The dust-corrected Millenium Galaxy Catalogue (Driver et al. 2007b) provides a field galaxy 
comparison for the Coma component luminosity functions. 

2.2. Structure and scaling laws of dwarf galaxies 

Early-type galaxies exhibit well-defined empirical correlations between global galaxy 
parameters, such as luminosity, radius, surface brightness, color, velocity dispersion, and 
line strength indices (Kormendy 1985; Faber & Jackson 1976; Bender et al. 1992; Guzman 
ct al. 1993; Graham & Guzman 2003). These scaling laws can tell us a great deal about the 
physical processes operating during galaxy formation, of which major contributors are star- 
formation feedback (Efstathiou 2000), tidal interactions (Due et al. 2004), and interactions 
with the hot intergalactic medium (Babul & Rees 1992; Roberts et al. 2004; TuUy et al. 
2002). 

Despite being the most numerous galaxy type in nearby clusters, dE/dSO galaxies 
are among the most poorly studied due to their low surface brightness (21< /Xg <25 V 
mag/arcsec^). The scaling laws for dwarf (My > —19) ellipticals and spheroidals are some- 
what controversial, it is unclear whether dwarfs form a family of galaxies distinct from giants 
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(Wirth & Gallagher 1984; Kormendy 1985; Caon et al. 1993), or a continuous sequence with 
them (Caldwell 1983; Graham & Guzman 2003). The ACS survey will allow sophisticated 
parameterisation of the surface brightness profiles of galaxies to \^ ~ 21 [My ~ —14); and 
measurement of basic structural parameters such as half-light radius and Sersic (1968) index 
to \^ ~ 23 [My ~ —12). Such limits are similar to those for measuring velocity dispersions 
and line strengths with large ground-based telescopes. Substantial samples from multi-object 
spectrographs on 4 metre class telescopes down to S ~ 19 already exist (see section [1]). 

The study of the scaling laws in various environments also provides the key observational 
reference needed to test specific predictions of current theoretical models of dwarf galaxy 
formation and evolution. For instance, the galaxy harassment model predicts that dE/dSOs in 
the highest-density cluster regions should have steeper light profiles, the fraction of nucleated 
dE/dSOs should be higher, the fraction of any remaining disk structure should be lower, and 
they should have higher metallicities than those located in lower-density cluster environments 
(Moore et al. 1998). 

2.3. The effect of environment upon morphological components 

The high density in the Coma cluster core makes it an ideal place to investigate the 
morphology- density relation, in which the average bulge-to-disk {B/D) luminosity ratio in- 
creases with galaxy density (Dressier 1980). This is understood to result from transformation 
of spirals into early types through processes such as harassment (Moore et al. 1996). An 
example of this transformation in progress may be the Coma cluster galaxy CMP 3629 
(Graham, Jerjen, & Guzman 2003). Furthermore, there is a strong radial gradient in HI 
content in the spirals, those in the cluster core being severely depleted (Bothun et al. 1984; 
Bravo- Alfaro et al. 2000). This, together with similar results in Virgo, supports the notion 
of gas stripping and the subsequent cessation of star formation (Conselice et al. 2003). From 
a quantitative morphological analysis using high-quality ground-based images, Gutierrez et 
al. (2004) find that disks in the Coma core are 30% smaller than in the field for a given bulge 
size. 

Bulge-disk decomposition in the presence of nuclear components depends critically on 
spatial resolution (Balcells et al. 2003). The high resolution of ACS allows identification of 
disks, unobscured primary and nuclear bars, classical bulges, resonance starburst rings, com- 
pact disks, disky bulges, and pseudo-bulges (e.g., Kormendy & Kennicutt 2004; Athanassoula 
2005). While observations in the rest-frame J-band (F814W) may miss some highly obscured 
morphological features, a comparison of the optical properties across field and dense clusters 
will set important constraints on how environment influences the morphology and structure 
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of galaxies. The radial dependence of bulge and disk morphologies can constrain the roles 
of mergers (Aguerri et al. 2001) and of disk truncation processes. 

2.3.1. Stellar Bars and Disks 

Stellar bars redistribute angular momentum in disk galaxies, driving their dynamical 
and secular evolution (Erwin et al. 2003; Kormendy & Kennicutt 2004; Jogee et al. 2004, 
2005; Erwin 2005). Many barred galaxies host molecular gas concentrations of up to several 
thousand Mq pc~^ in the inner kpc, intense circumnuclcar starburst rings, and disky bulges 
(e.g., Jogee 1999; Jogee et al. 2005). Bars are abundant in the field out to 2; ~ 1.0, corre- 
sponding to the last 8 Gyr (Elmegreen ct al. 2004; Jogcc ct al. 2004; Zheng et al. 2005), 
however wc know much less about barred disks in dense chistcrs. Varcla ct al. (2004) note 
that bars are twice as frequent in perturbed galaxies as in isolated galaxies, and Elmegreen 
et al. (1990) note a higher frequency of bars in binary galaxy systems compared with iso- 
lated galaxies, so there is some evidence for variation with environment. The closest to a 
systematic environmental survey is that of van den Bergh (2002), who found the frequency 
of bars in spiral galaxies in clusters and in the field to be similar, although he did not probe 
the richest cluster environments. 

Bars and similar structures (e.g., nuclear bars, rings, and spirals) are detectable in HST 
images if they have diameters ~ 3-5 FWHM (Sheth et al. 2003; Jogee et al. 2004; Lisker et 
al. 2006). At the resolution of the Coma survey they will be detected down to sizes of r ~ 150 
pc. This will enable identification and characterization of unobscured large-scale bars across 
the Hubble sequence (e.g., Erwin 2005), and across different environments. In moderately 
inclined galaxies, the survey can identify compact disks and disky "pseudobulges" (Erwin 
et al. 2003; Kormendy & Kennicutt 2004; Athanassoula 2005; Jogee et al. 2005). Coma is 
an excellent case study for the significance of inner/nuclear components in SO and spiral 
galaxies, and a key reference for comparison with studies of nuclear bar and ring frequencies 
in local field samples (e.g. Erwin & Sparke 2002; Laine et al. 2002; Knapen 2005). 

The surface brightness profiles of stellar disks fall into three classes (Erwin, Beckman, 
& Pohlen 2005; Pohlen & Trujillo 2006; Erwin, Pohlen, & Beckman 2007): classic single- 
exponential profiles, downward-bending "type II" profiles (Freeman 1970; "truncations"), 
and upward-bending "type III" profiles ( "antitruncations" ) . Type II profiles are common in 
early-type barred galaxies, and appear linked to the Outer Lindblad Resonance of bars; type 
III profiles, on the other hand, appear anti-correlated with bars, and may be a signature of 
minor mergers (Younger et al. 2007). A comparison of barred SO-Sb galaxies in the Virgo 
Cluster and the local field shows a dramatic dichotomy; type II profiles are common in the 
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field but essentially absent in Virgo (Erwin 2007, in prep). If this can be confirmed for the 
Coma Cluster, it points to a clear role for the cluster environment in shaping the outer disks 
of SO's and early-type spirals, and suggests an important test for distinguishing models of 
SO formation in clusters from the field. 



2.4. Colour gradients and stellar populations 

The steepness of the radial color gradients of dwarf and giant galaxies directly reflects 
their merging history: monolithic collapse imposes an initially steep negative metallicity 
gradient (Carlbcrg 1984), which will be progressively diluted by subsequent major mergers 
(Bekki & Shioya 1999; Kobayashi 2004). Semi-analytic models of hierarchical galaxy forma- 
tion predict differences in the merger history as a function of galaxy mass and environment 
(Cole et al. 2000). The Coma survey provides an unbiased sample of cluster dwarfs, so that 
the global scaling relations of different dwarf subtypes (dE, dSO, dE-N, dirr) can be rehably 
determined. Coma ellipticals have metallicities that seem to correlate with galaxy mass, 
which is broadly consistent with monohthic collapse models (e.g. Forbes et al. 2005). In 
contrast, dwarf ellipticals and spheroidals display a variety of color gradients, even positive 
in some places, which is an indication of recent star formation (van Zee et al. 2004). 

Complemented by X-band and IFU observations, the ACS images can be used to in- 
terpret the observed color distributions in terms of ages and metallicities, using photometric 
techniques (James et al. 2006) and line indices (e.g. Poggianti et al. 2001; Mehlert et al. 
2003; Sanchez-Blazquez et al. 2006). 

2.5. Globular clusters and UCDs 

The Coma cluster hosts several bright early- type galaxies, which are known to have 
very rich globular cluster (GC) systems (Ashman & Zepf 1998), and many less massive 
galaxies with their own star cluster systems that were assembled and evolved in the dense 
Coma cluster environment. There have been HST/WFPC2 studies of the GC systems of a 
small number of bright Coma ellipticals: IC 4051 (Baum et al. 1997; Woodworth & Harris 
2000), NGC 4881 (Baum et al. 1995) and the cD galaxy NGC 4874 (Kavelaars et al. 2000; 
Harris et al. 2000). Marin- Franch & Aparicio (2002) find a wide range of specific frequency 
Sn values among the brighter galaxies in Coma, using the Surface Brightness Fluctuation 
(SBF) technique on ground-based data. 

The ACS Virgo cluster survey has provided a number of important new results on the 
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GC populations of early- type galaxies in the sparser Virgo cluster (e.g. Jordan et al. 2005; 
Peng et al. 2006a; Mieske et al. 2006). At the distance of Coma it is only possible to 
resolve half-light radii down to ~20 pc (0".04). Moreover, the peak of the Globular Cluster 
Luminosity Function (GCLF) in Coma is suspected to be at K = 27.88 ± 0.12 (Kavelaars 
et al. 2000), some 1.3 magnitudes fainter than our lO-o" photometric completeness limit for 
point sources. Nevertheless, if the GCLF has a Gaussian form with a = 1.4 magnitudes some 
18% of the GCs are expected to be brighter than this hmit, allowing us to study the color 
distributions, color-magnitude diagrams and spatial distribution of GCs in a much wider 
range of galaxies than in the targeted Virgo cluster survey (from My — -23.3 to -15), to a 
larger physical radius from the galaxy centers (up to 100 Kpc). Stacking the GC systems 
of galaxies in bins will allow study of their properties as a function of host galaxy type, 
luminosity and environment, and thus of the relationship between GC systems and their 
hosts, and the processes which lead to the formation of the rich GC systems of the massive 
galaxies (Pipino et. al. 2007) 



2.5.1. The spatial distribution of clusters and the epoch of reionization 

The old, metal-poor GC population in giant ellipticals may provide the best available 
tracer of their dark halos, and thus may be useful in testing the predictions of ACDM models. 
Moore et al. (2006) and Bekki & Yahagi (2006) find that the final radial distibution of the 
old GCs depends upon the redshift of truncation of GC formation, which in turn might 
depend upon the epoch of reinonization. If the truncation is earlier, then the final radial GC 
distributions will be steeper and more compact (Bekki & Yahagi 2006). 

Bassino et al. (2006) present a study of the globular clusters around NGC 1399, the 
dominant galaxy in the Fornax cluster. They note a good agreement with a projected dark 
matter NFW (Navarro et al. 1996) density profile. Abadi et al. (2006) test the radial density 
profile of GCs around the somewhat isolated galaxies M31 and the Milky Way, comparing 
with the numerical simulations. They report that such luminous halos are similar in shape 
to their dark matter counterparts. On the other hand Merritt et al. (2006) in a careful 
analysis of a set of ACDM simulations, find that the halo density distributions are better fit 
by the much earlier model of Einasto (1965, 1972, 1974) than by NFW or a number of other 
alternative models. 
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2.5.2. UCDs and transition objects 

Ultra-compact dwarfs (UCDs), at —13.5 < My < —11.5, have been detected from 
the ground in the Fornax and Virgo Clusters (Hilker et al. 1999; Drinkwater et al. 2000; 
Phillipps et al. 2001; Jones et al. 2006). HST imaging of the bright Fornax UCDs shows that 
they have half-light radii of 10—20 pc, larger than both Local Group GCs and typical dwarf 
galaxy nuclei, but smaller than any previously known dwarf galaxies (De Propris et al. 2005). 
Similar but fainter objects have since been found (Ha§egan et al. 2005 - their dwarf galaxy 
globular cluster transition objects, DGTOs). These overlap in luminosity with, the brightest 
GCs, leading to the question of whether UCDs and large GCs are indeed fundamentally 
different objects. UCDs and DGTOs appear to be a feature of denser environments, so it 
is important to quantify the population of this new type of object in Coma and compare it 
with the less rich Virgo and Fornax. 

Proposed scenarios see UCDs (see Hilker 2006, for a review) as the tidally 'threshed' 
remnant nuclei of former nucleated dwarf elliptical galaxies in the cluster core (Bekki et al. 
2001), the similar remains of late-type spirals with nuclear star clusters (Moore et al. 1998), 
more massive, perhaps intra-cluster versions of ordinary GCs (Hilker et al. 2004), merged 
(super)star clusters (Fellhauer & Kroupa 2002), products of massive starbursts during the 
merger of two gas- rich galaxies (Fellhauer & Kroupa 2005), or even left-over primordial build- 
ing blocks of central galaxies in the dense galaxy cluster environment (Rakos & Schombert 
2004). 

The point source magnitude limit of the Coma survey is deep enough to observe compact 
objects down to My ^ —9 mag, so we can detect both UCDs and DGTOs down to GC 
luminosities. At the distance of Coma, the half-hght radii of the brighter examples will be 
0.02 to 0.04 arcsec, corresponding to ~ 10 — 20 pc (Drinkwater et al. 2003). These sizes are 
similar to those of extended star clusters in the Virgo Cluster which have been successfully 
differentiated from point sources by Peng et al. (2006b) using ACS. They can be reliably 
distinguished from background galaxies at magnitudes y ~ 22 to 25 mag which have half- 
light radii typically between 0.2 and 1.0 arcsec (e.g. Roche et al. 1997). 

2.6. The nature of "E-|-A" galaxies 

Spectroscopic surveys (Caldwell et al. 1993; Poggianti et al. 2004) have identified several 
galaxies with post-starburst spectra in both the core and the infall region of Coma. Sug- 
gestions as to what triggered these bursts range from equal-mass galaxy mergers (Barnes & 
Hernquist 1992) to high-speed impact with shock fronts of infalling substructures (Poggianti 
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et al. 2004) or generally with the dense ICM (Tran et al. 2005). The key to distinguishing 
between these processes is the spatial distribution of the intermediate-age populations. It 
would be concentrated to the center if major mergers were involved (Barnes & Hernquist 
1992) while it would be either an extended or off-ccntcrcd phenomenon if the galaxy had 
experienced an interaction with the ICM. The Coma survey includes a number of "E+A" 
or "k+a" galaxies at a range of luminosity, and the morphological features and color maps 
generated from the two ACS passbands will provide a valuable diagnostic of the physical 
origin of the "E+A" phenomenon. 

2.7. NucleEir Stair Clusters and Central Massive Objects 

Massive black holes at the centers of spheroidal stellar systems correlate with the large 
scale properties of the host spheroid. Popular correlations involve: (i) spheroid velocity dis- 
persion (Ferrarese & Merritt 2000; Gebhardt et al. 2000); (ii) stellar concentration (Graham 
et al. 2001; Graham & Driver 2007); (iii) luminosity (McLure & Dunlop 2002; Marconi & 
Hunt 2003; Graham 2007); and (iv) mass (Marconi & Hunt 2003; Haring & Rix 2004). 

The relationship between nuclear star clusters and central massive objects is of increasing 
interest, although the formation of the former remains poorly understood. They are observed 
in about 80% of intermediate-luminosity, early-type galaxies (e.g., Ferguson & Binggeh 1994; 
Graham & Guzman 2003; Cote et al. 2006; Jordan et al. 2007; Balcells et al. 2007). Graham 
& Guzman (2003) find that the luminosity of nuclear star clusters in dEs correlates strongly 
with the luminosity of the host spheroid. This luminosity trend has also been shown to 
exist in the bulges of lenticular and early- type spiral galaxies (Balcells et al. 2003). It 
has also been proposed that they follow the same Mnucieus^^spheroid relation as defined by 
supermassive black holes (Ferrarese et al. 2006b; Balcells, Graham & Peletier 2007), and the 
same Mnucicus^Sersic index relation (Graham & Driver 2007). 

Graham & Driver (2007) suggest that the driving physical relation may be with the 
central stellar density (prior to core depletion) rather than a global property of the host 
spheroid. Given the known trend between central stellar density and host spheroid luminosity 
(e.g., Graham & Guzman 2003; Merritt 2006a), the relations with global properties may be 
secondary in nature. It is therefore important to examine the connection between the mass 
of the nuclear star clusters and the central stellar density (projected, and deprojected, 
po) of the host spheroid. 

Study of the luminosity function, color-magnitude relation, relationship between nuclear 
and host galaxy color, and the radial distribution within the cluster (Ferguson & Sandage 
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1989; Lisker et al. 2007), and comparison with the lower density Fornax and Virgo envi- 
ronments, may help shed some light on potential formation mechanisms. Measurement of 
spatial offsets between nuclear clusters and the outer isophotes may reveal an oscillation of 
the nucleus about the center of the potential, which will be of larger amplitude in shallower 
potentials (Binggeli et al. 2000). Furthermore, apparant double nuclei may in some cases be 
a sign of an edge-on nuclear disk (Kormendy 1988; Tremaine 1995), from which we may be 
able to determine the presence of a central massive black hole (e.g. Debattista et al. 2006). 

Giant elliptical galaxies display cores which are partially depleted of stars. A possi- 
ble mechanism is the wrecking ball action of supcrmassivc black holes (SMBHs) from the 
progenitor galaxies as they sink to the center of the newly-formed galaxy during a merger 
(Begelman et al. 1980; Ebisuzaki et al. 1991, Merritt et al. 2007). Cores may also be en- 
larged when black holes are ejected from galaxy centers by the gravitational wave rocket 
effect following coalescence of a SMBH binary (Redmount & Rees 1989; Merritt et al. 2004; 
Gualandris & Merritt 2007). 

Using the core-Sersic model (Graham et al. 2003; Trujillo et al. 2004) it is possible 
to quantify the sizes (and mass deficits, M^ef) of these cores, and to predict each galaxy's 
central black hole mass, Mbh, using the measured Sersic index n (e.g., Graham & Driver 
2007). Given that Mdef/Mbh scales roughly linearly with N, where N is the number of 

major dry mergers, such measurements can be used to place constraints on the dry merger 
history of such cluster galaxies (Graham 2004; Ferrarese et al. 2006a; Merritt 2006b), and 
to constrain iV as a function of galaxy magnitude. 

3. Survey Design 
3.1. Survey area. 

The ACS Wide Field Camera (ACS/WFC) has a field of view of 11.3 square arcminutes. 
The camera contains two 4096 x 2048 pixel CCDs, with an inter-chip gap of some 3 arcsec- 
onds. The capabihties of ACS at the time of our pre-anomaly observations are descrived in 
some detail by Pavlovsky et al. (2006). 

The aim of the survey is to provide a large sample of galaxies for study in a high density 
environment, and at the same time a comparison sample in a lower density region of the 
cluster for the specific science goals of examining the effect of environment upon morphology, 
structural parameters and stellar populations. Because of the density of confirmed cluster 
galaxies in the core (e.g. CoUess & Dunn 1996) we have adopted a tiling strategy in the core 
region, tiling a region of approximately 18 x 21 arcminutes, using 42 ACS pointings in a 7 
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X 6 pattern, with some overlap between adjacent tiles. An HST orientation of 282 degrees 
was chosen for these tiles in order to maximise the time for which the obeservations could be 
scheduled in two gyro mode (Sembach et al. 2006). A tile of the central mosaic was moved to 
the southern edge of the mosaic area, away from the star HD112877 (V=7.17) which would 
have a negative impact on ACS observations near its position. 

In the outer part of the cluster the density of known members is less than one per ACS 
tile, and we decided to target known members which can be used to address some of our 
primary science goals. The best studied region outside the core is the infall region around 
NGC4839 (Neumann et al. 2001) where there have been a number of photometric surveys 
(Komiyama et al. 2002) and spectroscopic surveys for poststarburst galaxies (Caldwell et al. 
1993), and line strengths and velocity dispersions of dwarf spheroidals (Poggianti et al. 2001; 
Matkovic & Guzman 2005, 2008; Cody et al. 2007). 40 further ACS pointings were defined, 
each containing one, or in most cases more, cluster members from these spectroscopic surveys. 
In some cases the orientation was left free in the HST Phase 2 submission, in others a range 
was defined in order to ensure that more than one target galaxy was included. Table |6] lists 
the positions and orientations of our survey tiles. In this Table Column 1 gives the HST 
Visit number within the program; column 2 the field title. Columns 3 and 4 give the field 
center RA and Dec, and column 5 the field orientation of the observation, or the orientation 
specified for visits which have not yet been completed. Column 6 lists the spectroscopically 
confirmed cluster members within the field, with identifiers: GMP from Godwin et al. (1983); 
K from Komiyama et al. (2002). Column 7 gives the date of the observation, and column 8 
the number of dither positions which have been observed by the time of the January 2007 
ACS failure. At this time, 21 of our 82 fields have been fully observed, and further 4 have 
either two or three dither positions observed. The survey is thus 28% complete. 

Figure [2] shows in the left panel the positions and orientations of the survey tiles, and 
the locations of NICMOS parallel observations. For cases where a range of orientations are 
specified in the Phase 2 submission the orientation shown is the midpoint of the range, and 
those tiles with free orientation are shown at 90 °. The positions and final orientations of 
the 25 fully or partially observed fields are shown in the right hand panel of this figure. 
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Fig. 2. — Left: Positions of ACS tiles for the survey proposed, overlayed on a DSS image 
of the core and infall region of the Coma cluster. The small squares represent the fields 
of NICMOS parallel exposures. Right: The survey as at the ACS failure of January 2007. 
Tiles shown have some or all of the proposed observations complete by this date. In this panel 
requested orientations are replaced by observed orientations if they differ. 
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3.2. Choice of passbands, exposure times and dither parameters 

Color information is essential to some of the most important aims of the survey, so we 
require two passbands. For these passbands we choose F814W, as the passband which will 
give the deepest data for structural and luminosity function studies, and F475W, which is 
a compromise between color baseline and speed. These passbands are close to Cousins I 
and SDSS g' respectively. Transformations from these ACS/WFC passbands to standard 
Johnson/Cousins passbands are given by Sirianni et al. (2005). Total system throughput as 
a function of wavelength for all ACS/WFC filters, from the work of Sirianni et al. (2005), 
is given on the ACS project website [J and we plot in Figure [3] the throughput for our two 
passbands. 



http: / /adcam. pha.jhu.edu/instrunicnt/photometry/ 
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Fig. 3 — Total throughput curves of ACS CCDs plus filter for F475W passband (left) and 
F814W passband (right), after Sirianni et al. (2005) 
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We defined a dithering pattern with a large (3 arcsecond) dither across the inter-CCD 
gap on ACS/WFC, plus a sub-dither to remove hot pixels, resulting in four positions over 
the two orbits. The large move was made between orbits. Each passband was observed 
at each dither position. The large dither point spacing was 3.011 arcseconds at a pattern 
orientation of 85.28°, and the hot pixel sub-dither spacing was 0.2412 arcseconds at a pattern 
orientation of 22.3239°. POS-TARG equivalent positions are given in Table [2l 
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Table 2. POS-TARG equivalent of dither pattern used 



Dither Position POS-TARG 



1 

2 
3 
4 



0.0 0.0 

0.222 0.092 

0.247 3.001 

0.469 3.093 
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Exposure times were 350 seconds in F814W and 640 seconds in F475W at each dither 
position. F814W frames are somewhat deeper than F475W, but it was not possible to 
distribute the time more towards F475W, as this would have resulted in F814W frames too 
short to dump the ACS buffer, and thus a substantial increase in overheads. In some visits 
the F475W exposure time for the final dither position could be slightly longer, due to savings 
in the reacquisition overhead (as compared with the original acquisition), resulting in a total 
F475W integration time of between 2560 and 2680 seconds. 

In Table[3]we present magnitudes and surface brightnesses for signal-to-noise ratio (S/N) 
= 10, calculated with the ACS Exposure Time Calculator (ETC). In these calculations the 
ETC is set to assume standard background conditions, and an elliptical galaxy spectral 
energy distribution. The ETC gives two values for S/N for a point source: that in the 
"default" aperture, which for ACS/WFC is a circle of 0.2 arcsec radius; and an "optimal" 
S/N, assuming that a PSF fitting algorithm is used. In columns 3 and 4 of Table [3] limits 
are presented in each passband for each of these S/N values to reach 10. 

The real completeness and reliability limits, and contamination with spurious sources 
near these magnitude limits will depend upon real background conditions, to which the halos 
of large galaxies will contribute, and upon real HST guiding performance. For point sources 
we have estimated the limits by injecting artificial point sources into our ACS frames from 
one of our outer fields, and recovering them with SExtractor. In column 5 of Table [3] we 
present the magnitude at which 90% of artificial injected sources are recovered. This limit 
is, as expected, between the two calculated limits, but encouragingly close to the "optimal 
S/N" calculated limits. 

For galaxies the limits will depend upon surface brightness and structural parameters. 
These will be discussed further in Paper II in this series (Hammer et al. 2008). 

In column 6 of Table [3] we present a calculated limit in mag arcsec"^, obtained (using 
the ETC) for a S/N = 10 for a uniform surface brightness region of 1 arcsec^. 

All magnitudes and surface brightnesses have been converted from the Vega magnitude 
system used by the ETC to AB magnitudes, using the AB magnitude of Vega in the chosen 
filters (e.g. Sirianni et al. 2005). 
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Table 3. Columns 3 and 4 - predicted 10a limits for optimal extraction of point sources; 
column 5 - limit for 90% recovery of injected point sources; column 6 predicted surface 
brightness (1 sq arcsec) limits; all in the AB magnitude system. 



Passband 


Exposure 


ETC Limit (S/N=10) 


Measured 


SB (mag arcsec ^) 




(seconds) 


"Default" "Optimal" 


90% recovery 


S/N=10 in 1 arcsec^ 


F475W 


2560 


g' = 26.75 g' = 27.6 


g' = 27.55 


g' = 25.8 


F814W 


1400 


Ic = 25.95 Ic = 26.8 


Ic = 26.65 


Ic = 25.0 
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4. ACS Data Reduction 

The ACS data processing was performed at STScI. It involved a dedicated pipeline based 
on a wrapper script using PyRAF/STSDAS modules that performed CCD data reduction 
and cosmic ray cleaning, as well as the combination of individual dithered images using 
MultiDrizzle (Koekemoer et al. 2003), which makes use of the Drizzle software (Pruchter 
& Hook 2002) to remove geometric distortion and map the input exposures onto a rectified 
output frame. We describe details of the ACS data reduction below. 



4.1. Initial CCD Reduction 

As the observations of each HST visit were obtained, the data were run through the 
STScI on-the-fly-reprocessing (OTFR) pipehne for ACS data. This pipehne reduction in- 
volves basic CCD reduction by means of the IRAF/STSDAS task calacs, which performs 
bias subtraction, gain correction, dark subtraction, fiat fielding, and the identification of 
bad pixels. Finally, the OTFR pipeline identifies the exposures taken through each filter 
(F475W vs. F814W) and combines them using the OTFR version of MultiDrizzle for ACS 
data. These first-pass pipeline products were used for quicklook purposes. Final images were 
created by a second-pass calibration procedure which contained the steps described below. 



4.2. Updated Reference Files 

A few weeks after the each visit's data were observed, we downloaded more accurate 
calibration reference files from the HST archive, in particular the up-to-date superbias and 
superdark reference files which were created from bias and dark exposures obtained contem- 
poraneously with the science data (see, e.g., Lucas et al. 2006). These files were used for the 
final calibrations. 



4.3. Bias Level Offsets 

In virtually all ACS datasets of our Coma survey, we encountered significant bias level 
offsets between the four quadrants of the MultiDrizzle output images. These offsets, which 
ranged from a few tenths to several Data Numbers, are due to randomly varying differences 
between the bias level in the parallel overscan region of a given ACS/WFC CCD chip and 
the level in its active region (see, e.g., Sirianni et al. 2003; note that ACS/WFC contains four 
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CCDs that are read out by their own amphfiers). A script was developed to measure these 
bias offsets between the four quadrants using iterative statistics on a large number of pairs of 
areas located close to (and symmetrically relative to) the quadrant-to-quadrant borders. The 
measured offsets were subtracted from the appropriate quadrants of the superbias reference 
file being used in the final image reduction of each set of exposures with a given filter. 

4.4. Image Registration 

In many cases, the spatial registration of images taken during each HST visit of this 
Coma treasury program was less accurate than expected, even though all images were 
planned to be taken without any change of prime guide star, telescope position or roll angle. 
This problem was due to two compounding circumstances: (i) several visits were forced to 
use single-star guiding due to unforeseen problems with the HST Guide Star Catalog ver- 
sion 2 in the region of the Coma cluster, (ii) the observations of our program unearthed a 
software bug in the ground system of the two-gyro guiding mode of HST for guide star re- 
acquisitions when single guide stars were used. The result of these issues was that significant 
spatial offsets are present between images (and their world coordinate system (WCS) header 
keywords) taken in different HST orbits, even if taken during one and the same visit. Since 
it is crucial to align images to better than ~ 5 — 8 milliarcseconds (i.e., < 0.15 ACS/WFC 
pixel) in order to achieve accurate cosmic ray rejection among the separate exposures within 
a visit and to retain satisfactory point spread functions (PSFs) in the combined image, we 
registered all images within a visit to a common astrometric grid as follows. 

First, we ran MultiDrizzle on every set of associated images (i.e., images taken in one 
visit with the same filter) using the driz_separate = yes setting. This created so-called 
"singly drizzled" images, which are the individual input science images in the same format 
as the default output files from MultiDrizzle, i.e., after correction for geometric distortion 
using the WCS keywords in the image header. SExtractor (Bertin & Arnouts 1996) was 
then run with a signal-to-noise threshold of S/N = 10 on each single drizzled image. The 
resultant catalogs were trimmed on the basis of object size and shape parameters, thereby 
rejecting the vast majority of cosmic rays, CCD artifacts, and diffuse extended objects. 
Further catalog trimming was done by visual inspection of the 2-dimensional structure of the 
catalog sources, only retaining compact sources (non-saturated stars and compact galaxies). 
Typically, 3-15 "good" sources remained available for proper image alignment. Using these 
object catalogs, residual shifts (called 'delta shifts' in MultiDrizzle nomenclature) and 
rotations (Ax, A|/, A^^) between the individual singly drizzled images and a reference image 
were determined using IRAF tasks xyxymatch and geomap. The reference image was always 



-28- 



taken to be the image observed first in the visit. The resulting delta shifts and rotations 
were then fed to a second run of MultiDrizzle using driz_separate = yes to verify the 
result of the alignment. This often prompted the need for a second iteration of the alignment 
process. The formal uncertainties of the alignment process stayed below 0.1 ACS/WFC pixel 
(as reported by the geomap task). 

4.5. Cosmic Ray Rejection 

Once the alignments between the individual images were determined, cosmic ray rejec- 
tion was done nominally in two steps. The first step was performed by using MultiDrizzle 
and its cosmic ray rejection algorithm, which uses a process involving an image containing 
the median values of each pixel in the (geometrically corrected and aligned) input images as 
well as its derivative (in which the value of each pixel represents the largest gradient from 
the value of that pixel to those of its direct neighbors; this image is used to avoid clipping 
bright point sources) to simulate a "clean" version of the final output image. For a typical 
640-second F475W exposure, ~ 100, 000 — 320, 000 pixels were affected by cosmic ray hits, 

i.e., ~ 0.6 2% of the 4096^ pixels. Thus, it is exceedingly rare for a pixel to be affected 

by cosmic rays during all four exposures (namely <3 pixels out of 4096^). The number of 
pixels affected by 3 cosmic rays out of 4 exposures was however significantly higher (up to 
~120 pixels). In the latter cases, the median value was replaced by the value of the one 
valid pixel if the median value was larger than that of the valid pixel by a 5a threshold. This 
process yielded satisfactory results for the bulk of the tiles, except for the central strip of 
pixels in each tile which only had contributions from two exposures due to the gap between 
the two CCD chips of ACS/WFC. Furthermore, the MultiDrizzle output images of the 
tiles which only had two or three successful exposures also contained several residual pixels 
affected by cosmic rays. For the central strips in the tiles with four exposures and the full 
tiles with only two or three exposures, we therefore used a second step of cosmic ray rejection 
by means of the lacosmic routine (van Dokkum 2001) using parameter values chosen after 
careful testing. 

4.6. Construction of Final ACS Images 

Before the final run of MultiDrizzle on the input images to produce the final image 
tiles using the cosmic ray masks created as mentioned in the previous subsection, we cre- 
ated sky variance maps for each individual exposure for the purpose of deriving appropriate 
weight maps for the final image combination. These maps contain all components of noise 
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in the images, except for the Poisson noise associated with the sources on the image. These 
maps were constructed from the sky values determined by MuhiDrizzle for each individual 
exposure, the fiatfield reference file, the dark current reference file, and the read-out noise 
values as listed in the image headers for each individual image quadrant (or read-out am- 
plifier). These sky variance images were then inverted and used as weight maps for the 
associated exposures. 

The final run of MultiDrizzle was performed by shrinking the input pixels by a factor 
0.8 at the stage when the PSF is convolved by the input pixel scale (i.e., the f inal_pixf rac 
parameter). This factor was arrived at after some experimentation with different values, and 
was a good match to the degree of subsampling induced by the dither pattern we used. The 
final output images were produced in the default unrotated frame of the ACS/WFC CCDs 
rather than with North up so as to facilitate further analysis with PSF matching procedures. 

5. NICMOS parallel observations 

NICMOS parallel observations are made with the NIC3 camera of NICMOS in the J 
and H bands. The scientific motivations for these parallel observations, apart from that they 
are free, are the investigation of the near-IR luminosity function (LF), to compare with the 
SPITZER LF of Jenkins et al. (2007); a study of intra-cluster GCs, where the combined ACS 
(g,I) and NICMOS (J,H) data can be used to separate intermediate-age GCs from old halo 
globulars (e.g., Puzia et al. 2002); and a study of the near-IR fundamental plane for dwarf 
galaxies, where the relative insensitivity of the IR luminosity to dust extinction reduces the 
scatter in in the FP. 

The NIC3 camera has a field of view of 51.2 arcsc square, with the 3 arcsecond dither 
to accomodate the ACS inter-CCD gap, the uniformly exposed area is approximately 48 x 
51 arcsec per visit. The NIC3 field centre is some 8.5 arcminutes from the ACS field centre, 
at a position angle which depends upon the telescope orentation. The NICMOS fields are 
thus fairly random positions in the cluster, and as can be seen from Figure [2] only four of 
them overlap the area observed with ACS before the instrument failure. NICMOS parallel 
observations were taken on all pre-SM4 visits, of which 21 have the full four dither positions, 
these are summarised in Table |H 
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Table 4. NICMOS parallel observations - field centres and exposures 



Visit 


RA 


Dec 


Exposure per filter (s) 


1 


13:00:53.82 


28:13:10.9 


2560 


2 


13:00:38.21 


28:13:10.9 


2560 


3 


13:00:22.61 


28:13:11.0 


1280 


8 


13:00:50.20 


28:09:59.8 


2560 


9 


13:00:34.61 


28:09:59.9 


2560 


10 


13:00:19.02 


28:09:59.9 


2560 


12 


12:59:47.81 


28:09:59.8 


1920 


13 


12:59:32.21 


28:09:59.8 


1280 


14 


12:59:16.61 


28:09:59.9 


1920 


15 


13:00:46.60 


28:06:48.8 


2560 


16 


13:00:31.01 


28:06:48.9 


2560 


18 


12:59:59.80 


28:06:48.8 


2560 


19 


12:59:44.21 


28:06:49.9 


2560 


22 


13:00:43.11 


28:03:38.0 


2560 


23 


13:00:27.40 


28:03:38.0 


2560 


24 


13:00:11.91 


28:03:38.0 


2560 


25 


12:59:56.31 


28:03:38.0 


2560 


33 


12:59:37.00 


28:00:27.0 


2560 


45 


12:58:49.90 


27:33:21.2 


2560 


46 


12:58:40.56 


27:31:19.1 


2560 


55 


12:57:22.73 


27:38:58.6 


2560 


59 


12:58:31.11 


27:20:26.4 


2560 


63 


12:56:26.40 


27:21:58.4 


2560 


75 


12:58:54.33 


27:54:33.0 


2560 


78 


12:57:37.50 


27:30:20.6 


2560 
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In 2560 seconds exposure per filter the calculated lOcx magnitude limits for optimal 
extraction of point sources in the AB magnitude system are J = 25.9 and H = 25.4. 

6. Data and Educational products 

Data products will include image data as processed by the ACS pipeline described above, 
and also object catalogs with a variety of structural parameters as described in Paper II and 
subsequent papers in this series. 

The calibrated ACS images produced by STScI have been ingested into Astro- WISEl^(Valentijn 
et al 2006, Valentijn & Verdoes Kleijn 2006). Astro- WISE connects the distributed research 
groups for the data analysis and is used to create the source catalogs for the ACS images 
and to model the surface brightness distributions of galaxies. The results will be publicly 
available via the Astro- WISE web services. 

ACS data will be associated in Astro- WISE with all COMA Legacy Survey products 
including derived products and ancilliary data (e.g. Subaru, INT and UKIRT imaging and 
multi-wavelength data). Astro- WISE also connects the survey to the EURO-VO. 

Coma Legacy Survey data products will also be made available as part of the MAST 
Treasury archive at STScjf]. Expected timescales for the distribution of data products are 
given in Table [51 

An education and public outreach (EPO) program has been designed whose goal is to 
share the valuable legacy dataset and the associated research with the public. The deliver- 
ables include: 

1. HST ACS images of Coma that will appear at StarDate Online Astronomy Picture 
of the Week; in ViewSpace, which is seen in museums across the country; and in the 
revised StarDate/Universo Teacher Guide; 

2. An activity and DVD targeted at 9th- 12th grade students; and revised Internet versions 
of the StarDate/Universo Teacher Guide, which will highlight research on the Coma 
Cluster from HST data. 

3. Five Stardate radio programs on the Coma Cluster in English, Spanish and German for 



^http: / /www. astro- wise. org/ 

■^http: / / archive. stsci.cdu/hst / tall.html 
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distribution on monthly compact disk to over 500 radio stations in the USA, Mexico 
and Germany. 



The EPO program is a coordinated effort between NASA, the University of Texas at 
Austin, the McDonald Observatory Education and Outreach Office, and the Space Telescope 
Science Institute (STScI) 
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Table 5. HST Data Products and Data Release Schedule 



Data Product 


Description 


Release date 


Raw Data 


HST pipeline calibrated images 


Immediate 


Public catalog 


Positions, magnitudes, geometry, color. 


May 2008 




from S Extractor 




Processed Image Data 


Co-added, CR-cleaned, with best 


December 2008 




reference files 




Structural catalog 


Multi-component structural analysis 


December 2008 


External data 


Redshifts and ground-based colors 


May 2009 
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GMP3376,GMP3424 






27 


Coma4_6 


12:59:17.10 


27:55:21.0 


282 


GMP3486,GMP3510,GMP3511,GMP3522,GMP3645 
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Coma4_7 


12:59:01.50 


27:55:21.0 


282 


GMP3719,GMP3782,GMP3855 






29 


Coma5_l 


13:00:31.50 


27:52:10.0 


282 


GMP2716 






30 


Coma5_2 


13:00:15.90 
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282 


GMP2716,GMP2753,GMP2852,GMP2897,GMP2910, 
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282 
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GMP4206,GMP4330,GMP4381 
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GMP4192,GMP4215 
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GMP5284,GMP5320,GMP5364 
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12:57:04.22 


27:31:34.5 


299.1 


GMP5102 


19/Jan/2007 


4 


56 


Outskirts_14 
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57 


Outskirts_15 


12:58:02.08 


27:33:54.0 


250-260 


GMP4447,GMP4469,GMP4579 






58 


Outskirts_16 


12:58:34.59 


27:33:43.3 


220-225 


GMP4117,GMP4156,GMP4255 






59 


Outskirts_17 


12:58:31.15 


27:11:58.5 


270.05 


GMP4135,GMP4259,GMP4294 


16/Jan/2007 


4 


60 


Outskirts_18 


12:58:15.55 


27:05:15.1 


any 


GMP4383 






61 


Outskirts_19 


12:58:48.72 


28:00:52.8 


265-323 


GMP3969,GMP3973,GMP4003 






62 


Outskirts.20 


12:56:36.03 


26:54:17.8 


any 


GMP5361 






63 


Outskirts_21 


12:56:29.80 


27:13:32.8 


265.0 


GMP5335,GMP5359,GMP5365,GMP5472 


21/ Jan/2007 


4 


64 
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12:57:56.51 


27:02:16.4 
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GMP4582,GMP4596 






65 


Outskirts_23 


12:57:21.08 


27:00:04.2 


220-260 


GMP4961,GMP4980,GMP4947 






66 


Outskirts_24 


12:57:08.88 


27:05:47.8 


200-290 


GMP5032,GMP5097 






67 


Outskirts_25 


12:56:37.84 


27:02:49.3 


230-250 


GMP5283,GMP5293,GMP5294,GMP5395 






68 


Outskirts_26 


12:57:14.03 


27:15:14.7 


any 


GMP5012,GMP5019 






69 


Outskirts_27 


12:59:12.21 


27:36:56.9 


265 


GMP3585,GMP3598,GMP3696 






70 


Outskirts_28 


12:57:21.38 


27:37:12.1 


240-250 


GMP4888,GMP4905,GMP4967,GMP4987 






71 


Outskirts_29 


12:56:29.63 


27:12:27.2 


250-295 


GMP5345,GMP5362,GMP5434 






72 


Outskirts_30 


12:55:38.57 


27:44:39.9 


284-286 


GMP5850,GMP5912 






73 


Outskirts_31 


12:57:03.95 


27:45:06.4 


225-235 


GMP5096,GMP5100 






74 


Outskirts.32 


12:58:25.87 


26:54:23.2 


245-255 


GMP4183,GMP4340 






75 


Coma7_7 


12:58:47.75 


27:46:12.7 


280.0 


GMP3949,GMP3958,GMP4035,GMP4060 


17/ Jan/2007 


4 


76 


Outskirts.34 


12:58:19.13 


27:45:44.6 
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GMP4341,GMP4364 
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12:58:03.52 


27:40:57.8 


any 


GMP4502,GMP4557 






78 


Outskirts_36 


12:57:10.80 


27:24:18.0 


314.52 


GMP5038 


27/Nov/2006 


4 


79 


Outskirts_37 


12:57:50.52 


27:38:38.8 


any 


GMP4632 






80 
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12:56:04.00 


27:09:03.0 
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GMP5676 






81 
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12:59:50.50 


27:44:48.9 


230 


GMP3071,GMP3132,GMP3176,GMP3192,GMP3195, 
GMP3211 






82 
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12:56:14.60 


27:30:23.2 


any 


GMP5546 







